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r*| ■ Abstract 
Oh! 

This paper presents a diffusion experimental study for ternary lysozyme- 
Na2S04-water system, from moderate precipitant concentrations into 
the supersaturated region and provides a complete set of four diffu- 
sion coefficients. These data are important in order to provide accu- 
' rate models of protein diffusion with applications in growth of protein 

C**> • crystals for X-ray diffraction studies. All three-component mutual- 

diffusion experiments reported here were performed by Rayleigh in- 
| terferometry at pH= 4.5, T= 25° C and at a mean lysozyme con- 

""^5 ■ centration (average of top and bottom solution concentrations) of 0.6 

mM (8.6 mg/mL). Four experiments, with different combinations of 
protein and Na2S04 concentration differences, were performed at each 
pH * 1 of five mean Na2S04 concentrations (0.1, 0.25, 0.5, 0.65 and 0.8 M), 

Qh' for a total of 20 experiments. In addition, we have measured dynamic 

light-scattering diffusion coefficients of the ternary system lysozyme 
' chloride-Na2S04-water. 
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Motivation 



The diffusion of protein is one of the fundamental processes occurring in 
biological systems, and it is also an important step in the crystallization 
mechanism. Obtaining protein crystal of good structural quality is often 
the main issue for three dimensional atomic resolution structure studies of 
biomacromolecules. Crystallization is an intrinsically non-equilibrium pro- 
cess, and concentration gradients will occur around the crystal. The protein 
crystallizes, reducing its concentration at the moving face of the growing 
crystal and creating a protein gradient between the bulk solution and the 
crystal. This gradient in turn causes multicomponent diffusive transport of 
protein and precipitant. Diffusion in protein crystal growth inevitably occurs 
under conditions for which no species has an uniform concentration raising 
the issue of multicomponent diffusion. 

Crystallization experiments conducted under microgravity conditions have 
yielded protein crystals that provided diffraction data of significantly higher 
resolution than the best crystal of these proteins grown under normal con- 
ditions [lj. The difference between microgravity and normal gravity is the 
magnitude of the buoyancy forces. Commonly, this difference is assigned to 
some consequences of the reduction of buoyancy driven convection in the 
micrigravity conditions. With convection, the lysozyme concentration in the 
bulk solution is more uniform j2] and the probability for nucleation of para- 
sitic crystals is strongly reduced. 

The complete description of an n-solute system requires an n x n matrix of 
diffusion coefficients relating the flux of each solute component to the gradi- 
ents of all solute components [3]. The importance of other species on protein 
diffusion follows from the one-dimensional flux relations |3J: 



in which the cross-term diffusion coefficients (off-diagonal elements (Dtj) v i ^ 
j) can be positive or negative. In ternary systems (n = 2), our case, the one- 
dimensional flux relations could be written as: 
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where Ji and J 2 - the protein flux and respectively salt flux, (D u ) v and (D 2 2)v 
- the main-term diffusion coefficients relating to the flux of component to its 
own concentration gradient, and (-Di 2 )„ and (D2\) v - the cross-term diffu- 
sion coefficients relating the flux of each component to the gradient of the 
other. For some systems |U|3], a cross-terms {Dij) v can have considerably 
larger magnitude than the main-terms (Da) v , as our measurements show 
for the Lys-Na2S04- Water system described here and in agreement with the 
results about Lys-NaCl- Water system 6J. The index v from the diffusion 
coefficients shows that the experiment were done under the assumption the 
volume change on mixing and changes in concentrations across the diffu- 
sion boundary were small. Consequently, with a good approximation, the 
measured diffusion coefficients may be considered to be for the volume-fixed 
reference frame j7] defined by: 



where Vi is the partial molar volume of the ith species, and the subscript 
denotes the solvent. 

The importance of multicomponent diffusion has been recognized in the crys- 
tal growth community [HI E| and a crystal growth model has properly ac- 
counted for multicomponent diffusive transport in lysozyme chloride-NaCl- 
water system jHjnni- The experimental multicomponent diffusion coefficients 
are essential for accurate modeling of protein transport, especially in view 
of the very large cross-term coefficient (D2i) v reported here. Moreover, the 
concentration of supporting electrolyte dependence of all the diffusion coeffi- 
cients should be important for supersaturation region and also for its directly 
contribution to the protein flux. 

Experimental section 

All the experimental work was performed al Texas Christian University, in 
the Chemistry Department. 

Materials. All the materials, solution preparation procedures, apparatus 
and density measurement procedures are described in the work j^j. We used 
a hen egg-white lysozyme, recrystallized six times purchased from Seikagaku 
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America. 

The molecular mass of the lysozyme solute, Mi, was taken as 14307 g/mol, 
and this value ^1] was used to calculate all concentrations after correction 
for the moisture and chloride content. Buoyancy corrections were made with 
the commonly used lysozyme crystal density [T2"l IT51 ITlj of 1.305 g/cm 3 . 
The molecular mass of water, M Q , was taken as 18.015 g/cm 3 and the molec- 
ular mass of Na2SC>4, M 2 , was taken as 142.037 g/mol. 
Mallinckrodt reagent HC1 (~ 12 M) was diluted by half with pure water 
and distilled at the constant boiling composition. This resulting HC1 solu- 
tion (~ 6 M) was then diluted (pH 1.2) and used to adjust the pH of solution. 

Preparation of Solutions. All solutions were prepared by mass with ap- 
propriate buoyancy corrections. All weighings were performed with a Met- 
tler Toledo AT400 electrobalance. Since the as-received lysozyme powder 
was very hygroscopic, all manipulations in which water absorption might be 
critical were performed in a dry glove box. Stock solutions of lysozyme were 
made by adding as-received protein to a pre-weighted bottle that had con- 
tained dry box air, capping the bottle, and reweighing to get the weight and 
thus mass of lysozyme. Water was added to dissolve the lysozyme, and the 
solution was weighed. An accurate density measurement was made and used 
to obtain the molarity of the stock solution. 

For ternary experiments, precise masses of Na2SC>4 were added to flasks con- 
taining previously weighed quantities of lysozyme stock solutions. These 
solutions were mixed and diluted to within 10 cm 3 of the final volume. The 
pH was adjusted, and the solutions were diluted to their final mass. 

Measurements of pH. The pH measurements were made using a Corn- 
ing model 130 pH meter with an Orion model 8102 combination ROSS pH 
electrode. The meter was calibrated with standard pH 7 and pH 4 buffers 
and checked against a pH 5 standard buffer. 

Density Measurements. All density measurements were made with a 
Mettler-Paar DMA40 density meter, with an RS-232 output to a Apple 11+ . 
By time averaging the output, a precision of 0.00001 g/cm 3 or better could 
be achieved. The temperature of the vibrating tube in the density meter was 
controlled with water from a large well-regulated water bath whose temper- 
ature was 25.00+ 0.01 ° C. 
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Free-Diffusion Measurements. For binary Na2SC>4-water and ternary 
Lys-Na2SC>4- Water we performed measurements for free-diffusion using the 
high-precision Gosting diffusiometer EE CZ] operated in its Rayleigh in- 
terferometric optical mode. The procedure for measuring binary (D2) v and 
ternary diffusion coefficients (Dy)„ were described in detail in the work j^J. 
In order to measure the four diffusion coefficients of the system, the exper- 
iments must be performed with at least two different concentration differ- 
ences at each combination of mean concentration fSJ EH EJj. For ternary 
experiments, for each pair of mean concentrations, two measurements were 
performed with aci = and the two with ot\ = 0.8 (a< - the refractive index 
fraction due to the zth solute jHj). 

In order to make the data analysis of the free-diffusion experiments we used 
the Fick's second law: 



for two solutes. We made the assumption that the concentration differences 
of the solutes across the initial boundary are small enough and the diffusion 
coefficients are constant [201 • Also the volume changes on mixing were neg- 
ligible, thus all the measured diffusion coefficients are given relative to the 
volume-fixed frame of reference defined by equation (4). 

Dynamic Light-Scattering Diffusion Coefficients. We measured the 
dynamic light-scattering diffusion coefficients Ddls for samples from all ter- 
nary experiments and for 0.1, 0.25, 0.5, 0.65 and 0.8 M mean concentrations 
of Na 2 S0 4 . 

Measurements were made using a Protein Solution DynaPro-801 TC molec- 
ular sizing instrument with a fixed scattering angle of 90° and the procedure 
was described in the work [Bj. This apparatus allowed us to calculate also 
the eigenvalues Ai and X 2 of the matrix of diffusion coefficients, the Ddls 
predicted by Leaist's theory j2I] and to compare with ours and also with 
interferometric value (Dn)„. 

The viscosity for ternary lys-Na2SO,4-water system was also measured us- 
ing an Ostwald viscosimeter. 




i = 



1,2 



(5) 



5 



Results 



Ternary diffusion experiments were performed on the lysozyme chloride- 
Na 2 SC>4-water system at pH= 4.5 and T=25° C. To obtain the four ternary 
diffusion coefficients we performed four experiments at the same mean con- 
centrations but with different values of ACj for the solutes. There were two 
experiments with ACi = and AC 2 7^ and two with ACi 7^ and AC 2 = 
at each mean Na 2 SC>4 concentration of 0.1, 0.25, 0.5, 0.65 and 0.8 M. The 
interferometric data for the diffusion coefficients (Dn)„, (D 22 )t,, (Di2)o and 
{D 2 i) v are reported in the Tables 1,2,3. 

Partial molar volumes values, V\, V 2 and V Q , were calculated for each 
component using eqs A- 7 (q = 2) and 5 in [22] and reported in the Tables 
1,2,3. 

Values of mean density d and Hi = (dd/dCi)T, P ,Cj,j^i in the Table 1 were 
calculated using densities of all eight solutions from each experiment set. 
Densities were assumed to be linear in solute concentrations respecting the 
equation [U]: 

d = d + H 1 (C l -d l )+H 2 (C 2 -C 2 ) (6) 

where C\ and C 2 are the averages of the mean concentrations for all four 
experiments in a series. 

Dynamic light-scattering experiments provided the diffusion coefficients 
Ddls for each Na 2 SC>4 mean concentrations and also the eigenvalues A« (i = 
1,2). We could make an direct comparison to the interferometric values 
(Dn)„, the smallest eigenvalues of the matrix of diffusion coefficients Ai which 
are reported in the Table 1,2,3 and the Dp^s data. The values for T>dls are 
approximately 4% higher than (Dn)„, and approximately 5% higher than 
the smallest eigenvalue Ai of the matrix of diffusion coefficients. 
The viscosity values for ternary solution are reported also in the Tables 
1,2,3 at each mean value of salt concentrations. 

The diffusion coefficients dependence on C 2 mean values are shown in 
the Fig.l. 

From Fig. 1 we could see also the coefficient (Dn)„, due to the gradient of 
lysozyme, corrected from viscosity. The values for main-term (Du)„ are 13% 
smaller than the values for the same coefficient but using the lys-NaCl-water 
system, at the same concentration 0. The cross-term (D 21 )„ for the flux of 
Na 2 SC>4 caused by the gradient of lysozyme chloride increases as the Na 2 SC>4 
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concentration increases. The values for (D 21 )„ are 25% smaller than the val- 
ues for the same coefficient using the NaCl as salt 6J. At 0.8 M Na 2 S0 4 , this 
term becomes 14 times larger than the Na 2 S04 main-term diffusion coefficient 
(D 22 )„. The cross-term (Di 2 )„ for the flux of lysozyme caused by the gradient 
of Na 2 S04 is small in comparison with all the other diffusion coefficients as 
it was expected, and 35% smaller than (D 12 )„ in the case of lys-NaCl-water 
[U]. Fig. 1 shows us also that the ternary main-term (D 22 )„, for the flux of 
Na 2 S04 caused by the own gradient of concentration, is smaller with 1.5% 
than the binary coefficient (D 2 )^, over the entire composition range, and 40% 
smaller than the (D 22 )u for ternary lys-NaCl- water j^J, for the same range of 
concentration. 

Conclusions We reported the complete set of multicomponent diffusion co- 
efficients for ternary lys-Na 2 S04-water system at concentrations high enough 
to be relevant to crystallization studies. We also made a comparison with 
the reported diffusion data about lys-NaCl-water |B] in order to understand 
the influences of the type of salt in the diffusion process and, implicitly in 
the crystallization process. 
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Ci(mM) 

6 2 (M) 

d(gcm~ 3 ) 

H^ltfgmol- 1 

H 2 (10 3 gmol- 1 

Vi(cm 3 mol^ 1 ) 

V 2 (cm 3 mol~ 1 ) 

V (cm 3 mol~ l ) 

/Wio-Ws- 1 ) 

\ 1 (10' 9 m 2 s- 1 ) 

(DnUlO-Ws- 1 ) 

(D 12 ) v (10- 9 m 2 s- 1 ) 

(D 21 ) v (10- 9 m 2 s- 1 ) 

(D 22 ) v (10-'mls^) 

r)(cp) 



Table 1 

0.6000 
0.1000 
1.012189 

4.286 
0.12500 
10050 
17.090 
18.067 
0.11970 
0.11700 
0.1169± 0.0001 
-0.000013± 0.000001 
2.49± 0.01 
0.9661± 0.0001 
1.042395 



0.6000 
0.2500 
1.030671 

1.149 
0.12230 
10182 
19.780 
18.058 
0.11230 
0.10840 
0.1090± 0.0001 
0.000108± 0.000001 
4.14± 0.01 
0.8826± 0.0001 
1.110034 
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Table 2 



Ci(mJVl) 


0.6000 


0.6000 


C 2 (M) 

7/ — 3\ 

d[gcm °) 

TT { 1 rV^ 7 1 

zu hn3«™^/ — 1 

-H 2 ( iU gmol 


0.5000 


0.6500 


1.060538 


1.078032 


4.104 


4.049 


U.ll (do 


U.llolU 


Vi(cm 3 mol~ r ) 


10209 


10257 


V2(cm 3 mo/ _1 ) 


24.720 


25.930 


V (cm 3 mol~ v ) 


18.026 


18.013 


DnLsil^rrfs- 1 ) 


0.10102 


0.09430 


A^lO-Ws- 1 ) 


0.09566 


0.08775 


(AO^io-^s- 1 ) 


0.0969± 0.0001 


0.0894± 0.0001 


( J D 12 ),(10- 9 m 2 S - 1 ) 


0.000132± 0.000001 


0.000134± 0.000001 


(^o^io-Ws- 1 ) 


6.75± 0.01 


8.26 ± 0.01 


( J D 22 ),(10- 9 m 2 S - 1 ) 


0.7791 ± 0.0001 


0.7294 ± 0.0001 




1.236427 


1.322355 



Table 3 



Ci(mM) 




0.6000 


<5 2 (M) 

d(gcm~ 3 ) 

H^ltfgmol- 1 




0.8000 




1.095745 




3.954 


H 2 {Wgmol~ l 




0.11502 


Vi(cm 3 mol^ 1 ) 




10341 


V 2 (cm 3 mol~ 1 ) 




26.980 


V (cm 3 mol~ 1 ) 




17.993 


D DLS (10- g m 2 s- 


') 


0.08683 


\ 1 (lO~ 9 m 2 s- 1 ) 




0.08021 


(D 11 UlO- 9 m 2 s- 
(D 12 ) v (l0- 9 m 2 s- 
(D 21 ) v (10- g m 2 s- 


- l ) 


0.0822± 0.0001 


- 1 ) 


0.000130± 0.000001 


- 1 ) 


9.50± 0.01 


(L> 22 )„(10- 9 mV 




0.6900± 0.0001 


rj(cp) 




1.418206 
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